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Abstract. A non-Hermitian complex symmetric 2x2 matrix toy model is used 
to study projective Hilbert space structures in the vicinity of exceptional points 
(EPs). The bi-orthogonal eigenvectors of a diagonalizable matrix are Puiseux- 
cxpanded in terms of the root vectors at the EP. It is shown that the apparent 
contradiction between the two incompatible normalization conditions with finite 
and singular behavior in the EP-limit can be resolved by projectively extending the 
original Hilbert space. The complementary normalization conditions correspond 
then to two different affine charts of this enlarged projective Hilbert space. 
Geometric phase and phase jump behavior are analyzed and the usefulness of the 
phase rigidity as measure for the distance to EP configurations is demonstrated. 
Finally, EP-related aspects of VT— symmetrically extended Quantum Mechanics 
are discussed and a conjecture concerning the quantum brachistochrone problem 
is formulated. 
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1. Introduction 

A generic property of non-Hermitian operators is the possible occurrence of non-trivial 
Jordan-blocks in their spectral decomposition [1 . For an operator -ff (X) depending 
on a set of parameters X = (Xi, . . . , X m ) £ M. from a space M. this means that, in 
case of a single Jordan block, two or more spectral branches Ai(X), .... Afe(X) may 
coalesce (degenerate) at certain parameter hypersurfaces V C M. under simultaneous 
coalescence of the corresponding eigenvectors <£>i(X), . . . , $a,(X): Ai(X c ) = •■• = 
A fc (X c ), $i(X c ) = ... = $ fe (X c ) = 6 (X C ) for X c e V. Spectral points of this 
type are branch points of the spectral Riemann surface and are called exceptional 
points (EPs) pp. At the EPs the set of the originally k linearly independent 
eigenvectors $i(X), . . . , $fc(X) is replaced by the single eigenvector 9o(X c ) and k — 1 
associated vectors 8i(X c ), . . . , Oj._i(X c ) which form a Jordan chain. Together they 
span the so called k— dimensional algebraic eigenspace (or root space) 6>,(X C ) = 
span[9o(X c ), . . . , ®k_i(X. c )] [UI2] so that the total space dimension remains preserved. 
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The construction extends straight forwardly to the presence of several Jordan blocks 
for the same eigenvalue A(X C ). In general, the degeneration hypersurface V C M. 
consists of components V a of different codimension codim V a = a. Higher order Jordan 
blocks require a higher degree of parameter tuning (they have a higher codimension) 
and a correspondingly lower dimension of the component V a - Due to the different 
dimensions of its components V a the degeneration hypersurface V = 1J V a itself has 

a 

the structure of a stratified manifold [3J . 

EPs occur naturally in quantum scattering setups [H [5] when two or more 
resonance states coalesce and higher order poles of the S-matrix form. Within the 
Gamow state approach such S-matrix double poles have been considered in [6l [7j [H [9] , 
whereas in the Feshbach projection operator formalism (one of the basic approaches 
to analyze open quantum systems) they naturally occurred in studies of nuclei |10j . 
atoms plj [12] and quantum dots [Ml [14] . EP-related crossing and avoided crossing 
scenarios have been studied for bound states in the continuum [TTJ [321 [THl E] as well 
as for phase transitions [HI HH [20l [21] . In asymptotic analyses of quasi-stationary 
systems EPs show up as hidden crossings [22] . EP-related theoretically predicted 
level and width bifurcation properties have been experimentally verified in a series of 
microwave resonator cavity experiments. In [23] . the resonance trapping phenomenon 
(width bifurcation) [10] has directly been proven. The fourfold winding around an EP 
has been found experimentally [24] in full agreement with the theoretical prediction 
[19l [25] and related studies [TH [26]. In [27] two- level coalescences have been associated 
with chiral system behavior. The geometric phase at EPs has been discussed in 
O [261 EZl EE1 [291 [301 El [32] . 

EPs play also an important role in the recently considered VT— symmetrically 
extended quantum models [331 (Ml [35] . There they correspond to the phase transition 
points between physical sectors of exact VT— symmetry and unphysical sectors of 
spontaneously broken -pT-symmetry [36l l37| l38 | f39 j [40]. 

Other, non-quantum mechanical examples where EPs play an important role are 
the optics of bianisotropic crystals [4Tj . acoustic models [42] . many hydrodynamic 
setups where EPs have been studied within pseudo-spectral approaches [13] as well 
as a large number of mechanical models [44] where they are connected with regimes 
of critical stability [45]. Recent results on magnetohydrodynamic dynamo models 
indicate on a close connection between nonlinear polarity reversal mechanisms of 
magnetic fields and EPs [46] . 

For completeness we note that the perturbation theory for systems in the vicinity 
of EPs dates back to 1960 [UJ (see also [5]) and that it is closely related to 
singularity theory, catastrophe theory and versal deformations of Jordan structures 
[48] . Supersymmetric mappings between EP configurations have been recently 
considered in [49] [50] . 

A correct perturbative treatment of models in the vicinity of EPs has to be 
built on an expansion in terms of root vectors (eigenvectors and associated vectors 
6i(X c )) at the corresponding unperturbed eigenvalue A(X C ) (see e.g. El [44]). For 
X ^ V (away from the EP at X c and from other EPs) the operator -ff(X) has a 
diagonal spectral decomposition with corresponding eigenvectors $i(X). Choosing 
the normalization of these eigenvectors away from the EP and without regard to the 
expansion in terms of root vectors leads to a divergence of the normalization constants 
in the EP-limit X — * X c . The diagonalization break-down at X c , the occurrence of 
the Jordan block structure and the singular behavior of the eigenvector decomposition 
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are generic and were many times described in various contexts (see e.g. |51U52| ). The 
natural question connected with the fitting of the root- vector based normalization and 
the diagonalizable-configuration normalization (and related controversial discussions 
on their physical interpretation [511 I52| ) is whether and how the singular behavior 
affects the projective Hilbert space structure of quantum systems. 

In the present paper we answer this question by resolving the singularity with the 
help of embedding the original Hilbert space Ti = C" into its projective extension 
CP 71 instead of projecting it down to CP" -1 as in standard Hermitian Quantum 
Mechanics. Diagonal spectral decompositions and decompositions with nontrivial root 
spaces live then simply in different (and complementary) affine charts of CP n similar 
like monopole configurations of Hermitian systems have to be covered with two charts 
(North-pole chart and South-pole chart) of the unit sphere S 2 }53j . 

The basic construction is demonstrated on a simple complex symmetric (non- 
Hermitian) 2 x 2— matrix toy model. The consideration of complex symmetric matrices 
sets no restriction because by a similarity transformation any complex matrix can be 
brought to a complex symmetric form (see, e.g. |54[|55j V The Hilbert space notations 
for the 2 x 2— matrix model are fixed in section [21 In section [3l following [32J [44] 
we derive the leading-order perturbative expansion in the vicinity of an EP at X c 
in terms of root vectors and fit it then explicitly with expressions of the diagonal 
spectral decomposition at X ^ X c . Combining geometric phase techniques for non- 
Hermitian systems [28] with projective Hilbert space techniques from [56j we generalize 
the projective geometric phase techniques of Hermitian systems to paths around EPs 
(section HI) . The corresponding monodromy group is identified as parabolic Abelian 
subgroup of the orthogonal group 0(2, C) and evidence is given that vector norm 
scalings are only due to complex dynamical phases whereas geometrical phases are 
purely real-valued and norm preserving. In section [5] we consider an instantaneous 
(stationary type) picture of the system. Within such a picture, we resolve the singular 
normalization behavior by projectively embedding the Hilbert space H = C 2 CP 2 . 
We discuss the necessity for an affine multi-chart covering of CP 2 in order to 
accommodate diagonal-decomposition vectors and root vectors at EPs simultaneously. 
The usefulness of the phase rigidity as distance measure to EPs is discussed in section 
[6] In section [7] some EP-related aspects of VT— symmetric quantum models are 
discussed and a conjecture concerning the quantum brachistochronc problem [34[ 157] 
is formulated. The Conclusions [8] are followed by | Appendix A| where auxiliary results 
on Jordan structures of complex symmetric matrices are listed. 

2. Setup 

Subject of our consideration is the behavior of a quantum system near a level crossing 
point of two resonance states — supposing that for an N— level system the influence 
of the other N — 2 levels is sufficiently weak. Under this assumption the setup can be 
modeled by an effective complex symmetric (non-Hcrmitian) 2x2 matrix Hamiltonian 

»-(: :)• »= ht - w 

The effective energies e 12 G C and the effective channel coupling w£C will in general 
depend on underlying parameters X = (A 1; . . . , X^) € M. from a space M. 
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For nonvanishing coupling u> ^ the Hamiltonian can be rewritten as 

H = E <g> I 2 + u> ( Z x } z \ (2) 
with I2 denoting the 2x2 unit matrix and 

Eo:=\(e 1+ e 2 ), Z := ^ . (3) 

In this representation the eigenvalues E± and eigenvectors $± of H take the very 
simple form 



E ± = E ±cj\/Z 2 + 1 (4) 

and 

* ± -(_ z± ^ TT ) C± , c ± eC*:=C-{0} (5) 
which makes the branching behavior most transparenlPl. From the overlap 

($+!$_) = $; T $_ (6) 



1 + \Z\ 2 - \Z 2 + II + 2Im (Z*\/Z 2 + 1 



one reads off that ($+|$_) = holds only for ImZ = and that for general Z G C the 
two states $+ and $_ are not orthogonal ($+|$_) 7^ 0. Following standard techniques 
[4"3] for non-Hermitian operators, we consider a dual (left) basis S± bi-orthogonal to 
*± 

(1T+ - ^JHi = 0, (Hfcl^ocfa, fe,/ = ±. (7) 

For complex symmetric H it holds H± oc so that the most general ansatz for the 
right and left basis vectors $± and H± can be chosen as 

*±=c±x±, S±=d±x±7 c±,d±£C* (8) 
x± : = ( 7 ^ . /55-TT J • ( 9 ) 



-Z±v*Z 2 + 1 _ 

The bi-orthogonality 

<S±|<I> T ) =d ± c TX l XT =0 (10) 

is ensured by the structure of x± an d holds for any value of the parameter Z G C. 
A normalization (S±|$±) = 1 would set two constraints on the four free scaling 
parameters c± , d± G C* 

(S±|$±) = d t c±x^X± = l, (11) 

so that the system would still have two free parameters which should be fixed by 
additional assumptions. Subsequently, we will mainly work with the bi-orthogonality 
properties of the vectors <I>±, H± and fix their normalization only when explicitly 
required. 

Due to the arbitrary scaling parameters c±,d± G C* of the right and left 
eigenvectors $±, 3± G H = C 2 (HJ) it is natural to consider equivalence classes of such 
vectors defined by corresponding lines 7r($±), 7r(S-|-). These lines form the projective 



1 ) The fact that <E>± depends only on the single parameter Z reflects the property that after rescaling 
the energy by l/u> and shifting it by —Eo/u> (these transformations do not affect the eigenvectors) 
the Hamiltonian J2} depends only on Z. 
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Hilbert space F(H) =H*/C* = CP 1 3 7r($±),7r(3±) [SSIIMIESI, where H* := H-{0} 
denotes the original Hilbert space with the point at origin {0} = (0, 0) deleted to 
allow for a consistent definition of P(W). The space P(W) is covered by a single chart 
of homogeneous coordinates (zq,zi) t £ 7i and two complementary charts of affine 
coordinates t/o 9 (1, Zi/zo), z ^ and J/i 9 (zo/^i, 1), z\ ^ 0. Comparison with the 
structure of the auxiliary vectors \± ® shows that the x± can be straight forwardly 
re-interpreted as points of the projective space CP 1 described by the affine coordinate 
over Uq: x± ~ i"(^±)- The vectors {$±,£±1 themselves can be understood as 
sections of the natural line bundle L = {(p, v) £ F(H) x H\ v = cp 7 c £ C*} [S3], i.e. as 
$ ± = 7r($±)®c±, E± = 7r(S±)®c£|., where 7r denotes the projection tt : TL* —> ¥(H). 
The bundle structure is locally trivial tt" 1 ^) ps J7 x C* 3 $± [610 



3. Jordan structure 



At an EP, the two eigenvalues B± coalesce E + = E- = Eq, According to (U|), this 
happens for Z 2 = — 1 and Z = Z c := ±i and via JSj) it is connected with a coalescence 
of the corresponding lines 7r($+) = 7r($_) =: 7r($o) encoded in 

*+=*-=» : =(-U = (i)- <i2) 

This means that the eigenvalue E has algebraic multiplicity n a (E ) — 2 and geometric 
multiplicity n g (E n ) — 1 and by definition the level crossing point is an EP of the 
spectrum. The bi-orthogonality (TIT))) of &± and S T is compatible with the coalescence 
of the lines due to the vanishing bi-norm XoXo — 0, i.e. the isotrop\F)l of xo-, ~ 
a generic fact holding for the (geometric) eigenvector at any EP [2j [44]. We note 
that the coalescence 7r($+) = 7r($_) = 7r($o) still leaves the freedom for the vectors 
<J>+ = c+xo and $_ = c-Xo of being two different sections <!>+ ^ $_ of the same fiber 
7t($ ) x C* over 7r($ )- 

The right and left eigenvectors $o, ^-o at the EP are supplemented by 
corresponding associated vectors (algebraic eigenvectors) $i and Si defined by the 
Jordan chain relations [44] 

[H(Z c )-E Q I 2 ]<f> = 0, {H{Z c )-E I 2 }<S>i = <S>v (13) 
\H{Z C ) - E I 2 } + S = 0, [H(Z C ) - E I 2 } + S x - S . 

From the inhomogeneity of these Jordan chain relations it follows immediately that 
the root vectors 3>o and 3>i as well as S and Si scale simultaneously and in a linked 
way with the same single scale factor cq and d,Q , respectively. This is also visible from 
their explicit representation (|A.9[) derived in |Appcndix A| 

f 1 \ -x ( -Z r 



v -Z c J ' 1 y u v 1 

So = <rq*d* ( ~f c Ys 1 = aq*- 1 ^ ( _^ \ (14) 

pvf , — 

q := v2cj, Z c = ±i =: fii, co,do G C* . (15) 



V2 



2 ) For completeness we note that the (right) eigenvectors $± and the dual (left) ones H-j- could be 
understood as elements of a vector bundle ¥(7i) X F and its dual ¥(H) X F* with pairing in the fibres 
(.|.) : F* X F p > C (see, e.g. [58]). The details of this construction will be presented elsewhere. 

3 ) The vector xo behaves similar like a vector on the light cone in Minkowski space. 
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The simultaneous scaling means that the lines 7r(<j>o), ?r(So) at the EP should 
be interpreted as the one-dimensional components (projections) of two-dimensional 
planes which span the root spac^f)| &(Eq) 12 and which scale as a whole with a single 
scale factor. Such a higher-dimensional (complex) plane-structure goes clearly beyond 
the one-dimensional line structure of the projective space P(7i) (mathematically one 
should extend the natural line bundle of the original projective space to a more general 
projective flag bundle J62j[63]jf}| and underlines the fact that the state at an EP itself 
is not an element of the projective Hilbert space P(H) in its usual understanding. 

The basis sets {^o^i} and {So, Si} satisfy the well known bi-orthogonality 
conditions [44] 

<S o |$ o ) = (S 1 |$ 1 )=0 

<S o |$ 1 ) = (S 1 |$o)=doc o ^0. (16) 

Again, a normalization (So|<I > i) = (Si|$o) = 1 would only lead to a constraint 
do Co = 1 on the scale factors, but would not fix them completely. Due to this scaling 
freedom the single line 7r(<I>o) of a given Jordan structure, in general, still allows for 
different sections $o,a 7^ ®o,b of the corresponding fiber 7r(<i>o) x C* 9 $o,a, &o,b, 
7r($o,a) = 7r(*o,b) = 7r($o)- 

Let us now consider in detail what happens when the system approaches one of 
the critical values Z c = ±i. For this purpose we use the well-defined (but completely 
general and arbitrary) ansatz 

Z = Z c + e, |e|«l, eeC (17) 

and expand the eigenvalues J3J and the line defining vectors x± ® m terms of s. This 
gives the leading contributions to their Puiseux series representation [21 [44] in s 1 ^ 2 as 

E± = E ±e 1/2 AE + o(e 1/2 ), (18) 

AE := uj^/2Z~ c , 

(-k) ±eV, ( A) +o(e,/,) ' <19) 

Following [3H [33], we expand the eigenvectors <&±(Z) of the diagonal spectral 
decomposition in the same local e 1 ! 2 approximation in terms of the Jordan chain 
(root) vectors $o,i 

$ ± = $ + e 1/2 (6o*o + h*i) + o(e 1/2 ), (20) 

Z r 

b a =±-^AE, b 1= ±AE. (21) 



The coefficients &o,i are obtained by a two-step procedure. Substituting (fT7|) . (fT8j) . 
(|20p into the eigenvalue equation and making explicit use of the chain relations (fT3| 
yields b\ and leaves bo still undefined. The coefficient bo is found by comparing the 
line structure^)! of $± with x± m (HU)- 

It remains to match the fiber sections — what can be done in two ways. One 
may assume a single scaling coefficient Co of the root space given and consider the 



4 ) In the present simplest model ©(Bo) fills the whole Hilbert space H = C 2 . 

5 ) Indications that all the root vectors of a Jordan chain should scale simultaneously with a single 
scale factor were given, e.g., in [8] for Gamow vector setups with higher S-matrix poles. 

6 ) The term e 1 / 2 6o'I > additionally present in l|20| l in comparison with the corresponding result in 
1321 is due to the different choice of the root (Jordan chain) vectors f&o, The chain relation J 13D 
shows that the associated vector $1 is defined up to additional <t>o contributions and can be replaced 
by any linear combination $1 -f a^o, a £ C 
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coefficients c± of the sections $± as derived objects. This leads to the identification 
c + = c_ = <rqc . Apart from this option, one may assume the scaling coefficients 
c± as primary objects and given so that they may take different values c + ^ c_. 
Correspondingly the scaling factor Co of the root space will then be fitted to c± so 
that it will take two different values 

Co,± = c ± /(aq) . (22) 

Both constructions are possible and compatible with the smooth fitting of the line 
structure encoded in the EP-limiting behavior 7r($±) — * 7r($o). 

In a way similar to the above two-step procedure with subsequent fiber fitting the 
left eigenvectors can be obtained as 

5£ = S* + e l / 2 (b E* + 6 X S;) + 0(8^1 (23) 
d± = a*qZ c do^± . (24) 

Here, bo and b% are the same as in (|21|) and full compatibility with the bi-orthogonality 
conditions j2J as well as with (JSJ) is easily verified by direct calculation. In case 
of a single scaling factor do of the dual root space the coefficients d± will coincide 
d+ = d- = a*qZ c do- 

Combining ([2H)l and (|2"5|) one finds the limiting behavior of the inner products as 

(S±|*±) = 26 1 d ,±co,±e 1/2 + o(e 1/2 ) 

= ^ d±c±e i/ 2+o(£ i/ 2) . (25 ) 

Here, one has to distinguish two normalization schemes. If one assumes the root vector 
sets {<&□, ^i}, {So,£;l} normalized, e.g., with do,± c o.± = 1 or d c = 1 in (flU|) then 
the scalar product (S±|$±) of the eigenvectors in the diagonal spectral decomposition 
(see (|25[1 ) vanishes in the EP-limit. Starting, in contrast, from normalized eigenvector 
pairs {$±,S±} of diagonalizablc Hamiltonians as in (jTTJ) , i.e. with (S±|$±) = 1, 
then the scale factor products d±c± diverge as d±c± cx e -1 / 2 for e — ► 0. Both 
normalization schemes are possible and compatible with the smooth limiting behavior 
7r($±) — > 7t($ ) of the lines encoded in x±( E ~~ * 0) ~ * Xo [cf. IJHJ]. We see that 
this special behavior is only related to the fiber sections and not to the fibers (lines) 
themselves. The two incompatible normalization schemes simply indicate on the need 
for two complementary charts to cover the whole physical picture in the vicinity of 
a 2 x 2 Jordan block J^. One of these charts (we will call it the root vector chart) 
remains regular in the EP-limit, whereas the other (diagonal representation) chart 
becomes singular. 

The situation is similar to the two affine charts required to cover the Riemann 
sphere CP 1 . Starting from homogeneous coordinates (xq,Xi) £ CP 1 one has the two 
affine charts Uq 9 (l,Xi/xo), xq ^ and U\ 3 (xq/xi,1), X\ ^ 0. The mutually 
complementary affine coordinates z := Xi/xq € C 1 and w := xq/x\ G C 1 are then 
related by the well known fractional transformation w — 1/z so that the singular 
| z | — > oo limit in the z— chart corresponds simply to the regular w — > limit in the 
it;— chart. In other words, the two charts cover the North-pole region and the South- 
pole region of the Riemann sphere — a construction well known, e.g., from complex 
analysis and the description of magnetic monopoles |53j . 

Returning to the two-chart picture of the normalization we see that the original 
Hilbert space Ti = C 2 should be extended by the set of infinite vectors <E>± what can 
be naturally accomplished by embedding it into a larger projective space Ti. <—* CP 2 . 
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Correspondingly the fibers 7r($±) x C* should be extended as 7r(<I>±) x C* 
7r($±) x CP 1 . A detailed discussion of these structures will be presented in [M] . 
An explicit embedding construction for simplified setups with coinciding scale factors 
d± = c± is given in section [5] below. 



4. Geometric phase 

Following earlier studies [2JJ [28] [29] [30l [31] , geometric phases [65] of eigenvectors of 
non-Hermitian complex symmetric operators have been recently considered in [32] 
for paths in parameter space encircling an EP. The results showed full agreement 
with the phase considerations of [57]. In this section, we combine techniques for 
non-Hermitian systems |28[ 132] with explicit projective space parameterizations for 
Hermitian systems [56 to provide an explicit projective-space based derivation of the 
phase representation for non-Hermitian systems. Such an explicit reshaping of the 
results of [56j to non-Hermitian setups seems missing up to now. 

Following [28j [29[ [30] [32] we consider an auxiliary system with a general non- 
Hcrmitian Hamiltonian H{t) depending on a set of non-stationary parameters X(t) = 
[Xi(t),...,X m (i)] G M, H(t) = H[X(t)] and an EP hyper-surface V C M which 
is encircled by an appropriate loop T in parameter space M 3 V = {X(t), t € 
[0, T] : X(0) = X(T)}. The evolution of the quantum system is governed by a usual 
Schrodinger equation for the right eigenvectors 

»&*»(<) =H(t)* n (t), (26) 
and, due to the time invariance of the bi-orthogonal product 

{5 m (t)\$ n (t))=8 ma , (27) 
by a complementary evolution law for the left eigenvectors [28j 

id t Z m (t)=H+(t)Z m (t). (28) 

For an adiabatic motion cycle T C M. with Hamiltonian H[K(T)] = iJ[X(0)] the 
resulting eigenvector Q n (t = T) of _ff[X(T)] must lay on the same line as the initial 
& n (t = 0), i.e. it can only obtain an additional scaling factor which we parameterize 
as complex- valued phase 

$„(T) = e^ T ^ n (0), MT) G C. (29) 

Due to the preserved orthonormality (|27[) the corresponding left eigenvectors evolve 
as 

3 m (T) = e^( T )H ro (0). (30) 
The complex phase <j> n (T) can be split into a dynamical component [28] [56] 

e n{ T) = -( T (31) 



(S n (t)\$ n (t)) 

and the geometric phase 

7n(T) = <P n (T) - e n (T) . (32) 

Adapting the techniques of 56J wc calculate j n {T) in terms of explicit projective space 
coordinates. Setting 

$n(*) = c n (t) X n(t) =: [z (t), Zl (t)} T = z (t)[l,w(t)] T 

E m (t) = d* m (t) X * m (t) =■ [yo(t),yx(t)] T = y (t)[l,v(t)] T (33) 



Projective Hilbert space structures 



9 



(omitting in the projective space coordinates the mode indices to, n) one identifies 

<j> n {T) = -i\n[z {T)/z {Q)} (34) 
and obtains from ([26]) , ([31]) and ([32)) the differential 1-form of the geometric phase as 

.dz .y^dz Q + y'^dzx . v*dw 

"7= — « h« — t = l- . (35) 

z Vlzo + Vlzi l + v*w* 

Due to the symmetry ((8]) between left and right eigen-lines this simplifies to 



■ X T dx 



% 



d 1 = i^^ = ^d\YL(l + w 2 ). (36) 
XX 2 

Similar to results on Hermitian systems [5^ the differential 1-form ([551) i s independent 
of the coordinates zo and yo along the fibers and, hence, defines a horizontal connection 
over the projective Hilbert space of the system. The mere difference in the definitions 
of the projective structures is in CP 1 = S 3 /S 1 for Hermitian systems, whereas 
CP 1 = H*/C* for non-Hermitian ones [51]. 

Let us now apply the general technique to the concrete 2 x 2— matrix model |T|). 
Parameterizing the cycle around the EP by (TIT]) with 

£ = re la , a S [0, 2tt], < r < 1 (37) 

one reproduces the 1-forms of the geometric phases of [2"9"[ 

i [ i 

d'yj- — — dine — — da H — d\nr . (38) 
4 4 4 

In a similar way one obtains the same 1-forms for the corresponding left eigenvectors 
S-fc. Upon integration over a full cycle a(T) = a(0) + 2ir, r(T) = r(0) one finds 

7±(T)-7±(0) = ~. (39) 

The relation between geometric phases j± and the cycle phase a can be gained 
also directly from the structure of the sections <I>±. These sections may be arranged 
as columns of a diagonalizable 2 x 2— matrix 

:= [$+(«), $_(«)]. (40) 

The evolution along a cycle is then encoded in the transformation matrix W^a) = 
&(a) [3?(0)] which for small e with ^ \e\ <C 1 can be calculated from the 
representation (fT9|) as 

e^f 
2iZ c sin(f) e l f 

The elements W(a) form an Abelian parabolic subgroup P of the complex orthogonal 
group 0(2, C) D P (see, e.g., [62 EH [66]) 

W(a +/3) = W(a)W(/3) = W{f3)W{a) (42) 

corresponding to the mapping e %a S S 1 ~ l^(1) l_ > -P C 0(2, C). For full cycles 
a = 2nN, N G Z they yield the monodromy transformations [67) 

-i 



(41) 



VFo := W(0) = h, W x := W(2tt) = . ;/ 
VF 2 := W(4tt) = V^ 2 (2tt) = -I 2 , 

W 3 := W(6tt) = -W(2tt), W(8n) =I 2 = W . (43) 
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The geometric phase and the monodromy transformations show the typical 
four- fold covering of the mapping a i— ► 7 which was earlier described in [HO [551 123 122] 
and experimentally demonstrated in [24" 1. A cycle around the EP in parameter space 
has to be passed four times in order to produce one full 2tt— cycle in the geometric 
phase. A (non-oriented) eigen-line 7r(<£> =^ $0) S CP 2 is already recovered after 
two cycles ir(W2&) — 7r(— $) = 7r($) — similar to the eigenvalue E which for the 
2 x 2— matrix lives on a two-sheeted Riemann surface with the same two branch points 
Z c = ±i as the line bundle. For the isotropic limiting vector $0 at the EP it holds 
(due to CE2)) 

W(a)<S> = e~ M/4 $ (44) 

so that the parabolic subgroup P 3 W(a) can be identified as invariance group of the 
projective line at the EP 

7T (W(a)$ ) = 7T (e- M/4 $ ) = 7t($ ). (45) 

We note that despite the non-Hermiticity of the Hamiltonian H the geometric 
phase is purely real — as for Hermitian systems. Relations (|38|) . (|39|) show that 
possible imaginary phase contributions (which would result in a re-scaling of the 
eigenvectors $±) are cancelled by the closed-cycle condition r{T) = r(0). Hence, 
the non-preservation of the vector norm in non-Hermitian systems is induced solely 
by a complex dynamical phase e and requires the presence of the bi-orthogonal basis 
where a decaying behavior of the right eigenvector^ 

$„ cx e^"*"^*, ($„|$„) = ||$„|| 2 cxe- r " 4 (46) 

is necessarily connected with increasing vector norms of the dual left eigenvectors 

H m oc e-*-^*, ||S m || 2 cxe r ™' (47) 

so that indeed (S m |$„) = 5 mn . This behavior is well known from resonances and 
Gamow vector theory (see, e.g. [7])- 

Comparison of J5J) and (|46)l . ([47]) shows that the formal ansatz S m = for the 
eigenvectors of the complex symmetric Hamiltonian (cf. pTlfT3l[T4] ) can be used only 
for instantaneous eigenvectors at a single fixed t = to (which formally can be set to 
to = 0) as well as for the subclass of real symmetric matrices (when the system becomes 
Hermitian and norm-preservation of the eigenvectors holds). In contrast, for explicitly 
time dependent non-Hermitian setups it only holds S m (i) cx $„(t), i.e. the dual basis 
vectors necessarily live on complex conjugate lines (fibers) 7r[3 m (t)] = (7r[$ m (f)])* but 
with E m (t) ^ for t . 

Aspects of the parameter dependence of the phases and scalings in an 
instantaneous picture with S m = <&* n together with the explicit EP-limit e — > are 
subject of the next section. 



5. Instantaneous picture 

In modern quantum physics not only the properties of natural systems such as nuclei 
or atoms are of interest, but rather the design and functionality of artificial quantum- 
system-based devices plays an essential role. In many cases, for the understanding of 

7 ) For simplicity, we show the relations for stationary non-Hermitian Hamiltonians H with constant 
complex eigenvalues E n = e n + = const ; e n , T n £ R. 
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the dynamical features of these man-made quantum systems the time dependence is of 
minor interest. The properties of these systems are mainly governed by the position 
and number of EPs, i.e. the level crossing points in the complex plane, and their 
dependence on external control parameters. In this context it appears natural to study 
the parameter dependence of level energies and widths as well as the corresponding 
eigenvectors in terms of the instantaneous picture with S m = <&* n and c± — d± . This 
picture is compatible with the Hermitian limit when Imei^ = Imw = in ([1]) and the 
condition S TO = $*j is fulfilled by definition^. 

We have to distinguish the two possible normalization schemes — the root- 
vector based normalization (|16p with cIqCq = 1 or c?o,±co,± = 1 and the diagonal- 
representation based normalization (fTTj) with d±c±x T X — !• 

In the root-vector based normalization scheme the conditions c± = d± and 
^o,±co.± = 1 together with the two relations (|22|) and (|24|) imply (in leading-order 
approximation in e) co.± = g?o,± and, hence, co,± = K with k = ±1 (independently of 
the signs ± in the index of co,±) as well as c± — d± = nag. We see that in leading-order 
approximation in e the scaling factors c± = d± are rigidly fixed and independent of e. 
A geometric phase (necessarily induced via an e— dependence) appears incompatible 
with this normalization. 

Let us now turn to the diagonal-representation based normalization (jlip . In the 
EP-limit e — ► the normalization condition (fTTj) for the eigenvectors ([5]) yields 

2" 



1 = (H±|$±) =$T$ ± = 



i + utV^+i 



« T2Z c ^2ZTec 2 ± (48) 
and we find the expected divergent scaling factors as 

4 « T2- 3/2 Z- 3 / 2 e-^ 2 =► ci-e^ 1 / 4 . (49) 

On the one hand, (|49|) reproduces the local four-sheeted Riemann surface 



structure connected with the geometric phase (|38|) . (|4Tj) . i.e. a fourfold winding 
around the EP is needed to return to an eigenvector pointing into the same complex 
direction as a starting vector. (In contrast to the root-vector normalization scheme 
full compatibility with the geometric phase setup holds.) 
On the other hand, it leads to divergent vector norms 

||<i>±|| 2 = ($ ± |$±)«2|c ± | 2 «|2e|- 1 /2 (50) 

for e — > 0. As it was indicated in section [31 the corresponding singularity can be 
naturally resolved by embedding the original Hilbert space H ~ C 2 into its projective 
extension H CP 2 3 <fi = (uq, U\,U2) so that the set of infinite vectors becomes well 
defined. Interpreting the two components zq and z\ of the vector (fiber section) 

$ = c(l,w) = (z ,zi) e C 2 (51) 

as affine coordinates on the chart U 2 9 ( m , !LL , 1), u 2 ^ 0, U 2 C CP 2 

$ = (c, cw) ^> (c, cw, 1) (52) 



8 ) When compatibility with the Hermitian limit is not required, then the bi-orthonormalization 
constraints do.±co,± = 1 or d±c±x±X± = 1 fix only two of the four constants co.±,do.± or c±,d± 
and the remaining two can be chosen arbitrarily. For instance, one may set the eigenvector scaling 
factors as co,± = C ^ 1 or c± = 1 so that do,± = C -1 or d± = (x±X±) what would define 
instantaneous pictures not compatible with the Hermitian limit. 
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we can identify $ with the point <f> € CP 2 with homogeneous coordinates 

(f> = (u ,ui,u 2 ) = (l,w, c -1 ). (53) 

The singularity \c\ — > oo at the EP corresponds then simply to the point <fio — 
(1,10,0) S CP 2 with U2 = and we see that the affine chart U% € CP 2 is no longer 
appropriate for covering <fio. This is in contrast to the root-vector normalization scheme 
where c is fixed and the chart U2 remains suitable for the covering. Within the present 
diagonal-representation normalization, instead, 4>o should be parameterized in terms of 
affine coordinates corresponding to one of the charts^ Uq or U± with uq ^= or u± 7^ 0. 
Most natural for our representation (fBTj) . (fS"3")) is the affine chart Uq 3 (1, ^f") wruc h 
can be used for a suitable projective representation of the fibre sections $ 

$ » (1, c -1 ) = ( X T , c- 1 ) « (tt($), c- 1 ) . (54) 

Interpreting the normalization condition (|48|) as constraint on the affine 
coordinates of $ in the chart U2 

= $ T $ - 1 = 4 + 4 - 1 (55) 



in 



one immediately sees that it is equivalent to the conic (singular quadric 

u 2 a + u\ - u\ = (56) 

in homogeneous coordinates which cover the whole CP 2 . This conic remains regular 
at EPs which merely correspond to configurations with 1*2 = 0. In terms of 
(x T , c -1 )— notations it reads 

X T X-c- 2 -0. (57) 

It is clear that the conic construction is straight forwardly extendable to Hilbert space 
embeddings Ti — C™ CP™ of any dimension n. We arrive at the conclusion that 
the appropriate state space for open quantum systems in an instantaneous setting 
will be related to the projective extension CP n of the original Hilbert space Ti = C™ 
with states identified with conies X)fc=o u k ~ u n — 0- This is in contrast to Hermitian 
systems where it is sufficient to project the Hilbert space TL* — C™ — {0} down to 
the base space CP™ -1 , i.e. tt : H* -> F(H*) ~ CP" -1 . In non-Hcrmitian setups 
each fiber 7r(<i>) x C* should be supplemented by 00. This suggests to extend them 
to 7r(<E>) x CP 1 . From the above construction we see that the singular behavior with 
regard to the two affine charts is only related to the scale factors c € CP 1 , whereas 
7r($) behaves smoothly and regular. On its turn, this suggests to reconsider the model 
dependent physical interpretation of the eigenvector self-orthogonality (isotropy) and 
the corresponding diverging or non-diverging sensitivity in perturbation expansions 
like in [5TJ [S2] as result of divergent or non-divergent normalization constants. 

The Hilbert space extension Ti — C 2 e -» CP 2 together with the observed 
simultaneous scaling of the whole root space &\ obtained in section [3l the upper 
and lower triangular (parabolic subgroup type) structure of the &\— related matrices 
in (|A.7jl . (|A.8j) and the parabolic subgroup structure (|4l"1) at EPs provides one 
more indication that the natural structure at EPs is connected with projective flags 



9 ) A projective space CP™ 9 (zq, 21, . . . , z n ) is covered by n + 1 affine charts Uk 9 
(|2- ) . . , t f l j 5|±1, . . . , with z h 5^ (see, e.g., [59] [61]) in straight forward dimensional 

extension of the two-chart covering of the Riemann sphere CP 1 mentioned in section [3] 
10 ) For conies and quadrics in projective spaces see, e.g. 1591 1621 168[ . 
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[55] . A study of Jordan chain related flag bundles and the mappings between their 
complementary affine charts will be presented in [64] . 
Returning to the e — > limit in (|49p we see that 

4—i => ^-±< 

cl c_ 

=> $+ -> ±i$_ , (58) 

i.e. the two eigenvectors (fiber sections) $ + , <f>_ are phase-shifted one relative 
to the other by ±i when tending to their common coalescence line at e = 0: 
7r($ + ) = 7r($_) = 7r($o)- We note that this relative ±i phase-shift of the vectors 
$ + , is generic for models in their instantaneous picture and with d± — c± and 
normalization (S±|$±) = 1. 

A further result which immediately follows from (|4T)]) is the typical distance- 
dependent phase jump behavior in the vicinity of the EP. In a sufficiently close 
vicinity of an EP (|e| -C 1) any sufficiently smooth trajectory in an underlying 
parameter space can be roughly approximated by a straight line segment with an 
effective parametrization of the type 

e = e iao (p + is), s € [-so, s ] C M (59) 

where «o =const fixes the direction orthogonal to the effective trajectory in the 
complex e— plane and p is the minimal distance p = \e(s = 0)| of this trajectory 
to the EP. The parameter along the path is s G [— s ,so] C K. This parametrization 
gives: 

[e( S )]- 1/4 =e- i ^- ie W| £ (s)|- 1 / 4 
\e(s)\ = (p 2 + S 2 ) 1/2 

6(s) = iarctan(s/p) e (-7r/8,7r/8) (60) 

and we observe that the minimal distance p between the parameter trajectory and the 
EP defines the smoothness of the phase changes. The closer the path approaches the 
EP the more it will take the form of a Heaviside step function with jump height 7r/4: 

6(s) - . (61) 



0(s; P ^0)^ \ 



The phase jump behavior can be used as implicit indicator of a possible close location 
of an EP — a fact especially useful in numerical studies of systems with complicated 
parameter dependence, but where phases of eigenvectors can be easily extracted. 
Jumps ±7r/4 of wave function phases have been observed numerically in [5j5] for the 
model Hamiltonian ([T]) and in [T3] for the special case of a small quantum billiard. 
According to these results, the phases of the components change smoothly (as a 
function of a certain control parameter) in approaching the EP and jump by 7r/4 
at the smallest distance from this point. Other phase jump values are possible, but 
require especially tuned paths. 



6. Phase rigidity 

In numerical studies of man-made open quantum systems depending in a complicated 
way on several parameters X = (Xi, . . . ,X m ) G M. it is usually important to know 
how close a given configuration is located to an EP. EPs dominate the system behavior 
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also in their vicinities, spectral bands may merge at EPs [51] or the transmission 
properties of quantum dots (QDs) may become optimal at EPs [70]. A measure for 
the distance between a given point in parameter space and a closely located EP would 
provide a convenient tool for adjusting and tuning parameters so that a system may 
be 'moved' in parameter space toward to or away from this EP. 

In [70] it has been shown numerically that within the instantaneous picture 
($ = H*) an appropriate measure for the detection of EP vicinities is the fraction 

r = — — . (62) 

We note that originally similar fractions have been introduced in [7T] to describe the 
transitions between Hamiltonian ensembles with orthogonal and unitary symmetry 
in Hermitian quantum chaotic systems. There the square modulus |r| 2 was dubbed 
"phase rigidity". In our considerations of non-Hermitian systems we use this term in 
loose analogy for r itself. 

Decomposing $ into real and complex components $ = $ r + i$>i we find from 
the normalization that 

$ T $ = 1 = $J$ r - $T , $^$, t = (63) 
an hence that the norm is bounded below 

||$|| 2 = ($|$) = $^$ r + $f = 2$f $i + 1 > 1 . (64) 
The phase rigidity can be expressed as 

r = W G[0 ' 1] (65) 

where according to ([50]) for the EP-limit e — > holds 

rral&l^-fO. (66) 

The opposite limit r — > 1 is reached when the channel coupling to in the Hamiltonian 
(TTJ) vanishes, i.e. when the interaction between the two decaying resonance states 
tends to zero and any eigenvector can be taken purely real-valued in the instantaneous 
picture. 

Finally, we note that for certain quantum dot systems the phase rigidity r is 
closely related to the transmission properties of these systems. The capability of 
corresponding numerical studies (including the visualizations of transmission and 
phase rigidity 'landscapes' over parameter space) has been recently demonstrated in 



7. VT— symmetric models 

Toy model Hamiltonians of 2 x 2— matrix type have been often used as test ground in 
PT-symmetrically extended Quantum Mechanics (PTSQM) [33J [M] |35] . They can 
be obtained from non-Hermitian complex symmetric 2 x 2— matrix Hamiltonians by 

xl ) In equation (163 ft it can be set <t>jr<t> r =: cosh 2 and Qj^i ='■ sinh 2 ji. This hyperbolic structure 
shows analogies with the mass shell condition E' 2 — p 2 = m 2 of special relativity. The EP-limit 
<l?^ < I , r , &f&i — » co corresponds, e.g., to the light-cone limit where the vectors become isotropic — a 
fact which seems to play an important role in connection with the conjectured Hilbert space worm 
holes |34j related to the brachistochrone problem of VT— symmetric Quantum Mechanics (PTSQM). 
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imposing a VT— symmetry constraint. In a suitable parametrization they have the 
form 

H=( r f re -,e )> r,s,8eR (67) 

and commute with the operator VT 

[PT,H}=0, -P=(\ (68) 

Here, V is the parity reflection operator and T — the time inversion (acting as complex 
conjugation). The eigenvalues of H are 



E± = r cos(0) ± yjs 2 - r 2 sin 2 (0) (69) 
and the corresponding eigenvectors can be represented as [35] 

e ia/2 / 1 \ 



where 



V2 


cos(a) 




-ia/2 



y/2 cos(a) 



— e 



1^-) = -, x ( ia ]=:c-X- (70) 



r 



sin(a) = - sin(0) . (71) 



s 



With regard to the indefinite (Krein space type [35]) VT inner product (u, v) — VTu-v 
the vectors are normalized as 

(E±,E±) = ±1, (E ± ,E T )=0. (72) 

The indefinite VT inner product is then mapped by the dynamical operator C with 
[C, H]=0 and 

C = ^(^) ,} \ (73) 
cos(a) V 1 -*sin(a) / 

(see, e.g., (35J) into the positive definite (Hilbert space type) CVT inner product 
((u, u)) = CVTu ■ v which yields 

((E±,E±)) = l, ((E ± ,E T )) = 0. (74) 

Let us now reshape the model in terms of the EP-relevant notations of section [2] 
A simple comparison of |T]), ([3J with ([67]) . (|7ip shows that 

Z = ism{a) (75) 

and, hence, that 

cos(a) \ 1 ) 

and that the model is actually one-parametric with essential parameter Z . Together 
with ^ the parametrization (|7rj|) leads to a representation of the Hamiltonian (|67[) as 

F = + s cos(a)C , E = r cos(0) (77) 

and [C,H] = 0is fulfilled trivially. 

The compatibility of the VT and the CVT inner products f72"]). ([71)1 with the bi- 
orthogonality relations is ensured by the fact that for an eigenvector $ = c(l, 6) T 
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exact VT symmetry requires VT& cx <!> and, hence, c*b*(l, l/b*) T cx c(l, 6) T so that 
\b\ 2 — 1. For such vectors it holds VT& cx S* and due to the dynamically tuned 
C also CVT$> cx S*. As result one finds CVT^ k ■ $ ; cx VT<S> k ■ $ ; cx and 
full compatibility of the bi-orthogonality with the VT and CVT inner products is 
established. 

From ([75]) we see that possible EPs are solely defined by the value of a. From 
Z c = ±i we find the corresponding critical a c as 

a c = ±7r/2 + 2iV7r, NeZ. (78) 

Furthermore, it follows from (|7ip that a purely Hermitian model with 8 — nit, 
n G Z corresponds to a = iV7r, JV e Z. Exact VT— symmetry is preserved for 
a £ R - {tt/2 + ttZ}, and the corresponding models are parameterized by elements Z 
belonging to the purely imaginary straight line segment connecting the two EPs, i.e. 
by Z G ReZ = 0. 

According to ([70|) . at the EPs the eigenvectors lay on the same line tt(\E+)) = 
tt(\E-)) f» xo = (1, Z C ) T and their norms diverge for a — > a c like 

lll^±>ll 2 = <^±|^±> « ■ — V^r — °°- ( 79 ) 

| cos(a)| 



The operator C in (|73p shows the same singular behavior, i.e. the C— induced mapping 
between the Krein space and the Hilbert space breaks down at the EPs. In analogy to 
the singularity resolution presented in section [5] we may map the vectors \E±) G C 2 
into elements from the affine chart U2 C CP 2 corresponding to points e± G CP 2 with 
homogeneous coordinates 

\E ± )^e ± = (xl,c ± 1 ) . (80) 

The original normalization via VT inner product VT\E±) ■ \E±) = 1 acts then as 
generalized conic 

VT X ± ■ x± - (Tc^) Ci 1 = (81) 

which remains regular in the EP-limit a — > a c , but shows the typical EP-related self- 
orthogonality (isotropy) of the lines VT\± 'X± ~* 0- Again we arrive at the conclusion 
that the original Hilbert space TL — C 2 should be projectively embedded into CP 2 in 
order to accommodate EP-related singularities. 

Finally, we note that the recently uncovered solutions of the VT— symmetric 
brachistochrone problem with vanishing optimal passage tim e 134] occurs for a = 7r/2 
what according to (f78|) can be identified as an EP-regima 12 l This fact appears 
compatible with the results of 57J where a vanishing passage time was reported for 
arbitrary non-Hermitian Hamiltonians. In this regard it is natural to conjecture that a 
vanishing optimal passage time might be a generic EP-related feature of non-Hermitian 
systems not necessarily restricted to PTSQM models. 

8. Conclusion 

In the present paper we considered projective Hilbert space structures in the 
vicinity of EPs. Starting from a leading-order Puiseux-expansion of the bi- 
orthogonal eigenvectors of a non-Hermitian (complex symmetric) diagonalizable 
2 x 2— matrix Hamiltonian in terms of root vectors (algebraic eigenvectors) at an EP the 



J ) The corresponding state vector alinement without link to EPs was observed also in [72) . 
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normalization divergency of the eigenvectors in the EP-limit has been parameterized. 
It has been shown that the natural projective line structure related to the eigenvectors 
of the diagonal Hamiltonian has to be replaced at an EP by a higher dimensional 
projective structure in which all the root vectors of a Jordan block scale simultaneously 
with the same single factor. For a simplified setup with left eigenvectors equated to 
their complex conjugate right counterparts, the normalization divergency has been 
resolved by embedding the original Hilbert space 7i — C 2 into its projective extension 
H CP 2 . Eigenvectors normalized according to the diagonalizable Hamiltonian and 
eigenvectors with a normalization inherited from the root vector normalization live 
then merely in different (complementary) affine charts of CP 2 . The states themselves 
can be interpreted as conies in CP 2 . The line structure of the states behaves smoothly 
and independently of these charts and their possibly singular transition functions. This 
indicates on the possibility of a technically efficient description of the global behavior 
of the non-Hermitian system by factoring the eigenvectors in globally smoothly vary ing 
non-singular projective line components 7r(<&/-) and possibly diverging scale factora 13 ) 



With the help of the Puiseux expanded eigenvectors it has been shown that the 
geometric phase obtained on circles around EPs of complex symmetric Hamiltonians 
is purely real-valued and that the corresponding monodromy transformations are 
induced by an Abelian parabolic subgroup of 0(2, C). Furthermore the Puiseux 
expansion has been used to explain phase jumps which in prior work had been 
numerically observed in the vicinity of EPs. An analytical foundation for the usefulness 
of the phase rigidity as distance measure to EPs has been provided. Finally, a 
VT— symmetric model has been studied. It has been shown that the EP-related 
singularities show up not only in the normalization conditions of the eigenvectors 
but also in the dynamical symmetry operator C. The normalization singularity has 
been resolved via a projective extension of the original Hilbert space. From the 
singularity structure it has been conjectured that the zero passage time effect in 
the brachistochrone problem of non-Hermitian Hamiltonians might be a generic EP- 
related artifact. 
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Appendix A. Jordan normal forms for complex symmetric 2x2 matrices 

At the EPs with Z c = ±i =: \ii the matrix 



13 ) The question concerning the physical interpretation of diverging or non-diverging normalizations 
and the corresponding diverging or non-diverging sensitivity in perturbation expansions is highly 
model dependent (see e.g. 3 i 52_) and still requires a detailed investigation. 




M = PRJ 2 {0)R- 1 P- 1 . 



(A.2) 
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From the symmetry properties 

M = M T , J 2 (0) = S 2 J2(0)S 2 (A.3) 

with S2 = I ? n ) an d -P 2 = &2 one finds 



1 



and 



iJ= ( g -i )- g:=V^. (A.5) 

The elementary Jordan block J2(0) has right and left root vectors 80, 61 and ^0, ^1 
satisfying 

J 2 (0)6 = 0, J a (0)9i = 6 

J 2 T (0)* = 0, J 2 T (Q)* 1 = * - (A.6) 

The explicit solutions of these Jordan chains can be arranged as Toeplitz and Hankel 
matrices 



c o ci ^ iT , _ riTr lTl 1 ( d* 



and 



= [e o ,e l] = ( - - i, * = [* ,*i] = l £ ,S ) (A.7) 



*:=*&=( 5 £ ) • (A- 



d\ d* 



From the simplest realization of the bi-orthonormality condition \& + @ = S%, = 
I2 one finds the parameters c\ = d% = 0, d^c® = 1. Via similarity transformations 
©o,i ► ^0,1 = PRQos and $0,1 ► 3o,i = P ^0,1 one arrives at the root 

vectors of M 

$0 = crgco f V $1 = aq~ x CQ ( ^ c 
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